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ABSTRACT 

The  onset  times  and  peak  intensities  of  solar  energetic  particle  (SEP)  events  at  Earth  have  long  been  thought  to  be 
influenced  by  the  open  magnetic  fields  of  coronal  holes  (CHs).  The  original  idea  was  that  a  CH  lying  between  the 
solar  SEP  source  region  and  the  magnetic  footpoint  of  the  1  AU  observer  would  result  in  a  delay  in  onset  and/or 
a  decrease  in  the  peak  intensity  of  that  SEP  event.  Recently,  Gopalswamy  et  al.  showed  that  CHs  near  coronal 
mass  ejection  (CME)  source  regions  can  deflect  fast  CMEs  from  their  expected  trajectories  in  space,  explaining 
the  appearance  of  driverless  shocks  at  1  AU  from  CMEs  ejected  near  solar  central  meridian  (CM).  This  suggests 
that  SEP  events  originating  in  CME-driven  shocks  may  show  variations  attributable  to  CH  deflections  of  the  CME 
trajectories.  Here,  we  use  a  CH  magnetic  force  parameter  to  examine  possible  effects  of  CHs  on  the  timing  and 
intensities  of  41  observed  gradual  E  ~  20  MeV  SEP  events  with  CME  source  regions  within  20°  of  CM.  We  find 
no  systematic  CH  effects  on  SEP  event  intensity  profiles.  Furthermore,  we  find  no  correlation  between  the  CME 
leading-edge  measured  position  angles  and  SEP  event  properties,  suggesting  that  the  widths  of  CME-driven  shock 
sources  of  the  SEPs  are  much  larger  than  the  CMEs.  Independently  of  the  SEP  event  properties,  we  do  find  evidence 
for  significant  CME  deflections  by  CH  fields  in  these  events. 

Key  words:  acceleration  of  particles  -  Sun:  coronal  mass  ejections  (CMEs)  -  Sun:  flares  -  Sun:  particle  emission 
Online-only  material:  color  figures 


1.  INTRODUCTION 
1.1.  SEP  Events  and  Coronal  Holes 

Despite  their  important  space  weather  effects  (Bothmer  & 
Zhukov  2007;  Guetersloh  &  Zapp  2010),  the  occurrence  and 
timing  of  gradual  solar  energetic  ( E  >  10  MeV)  particle 
(SEP)  events  are  difficult  to  predict  (Crosby  2007;  Kahler  et  al. 
2007).  The  consensus  view  is  that  the  SEPs  are  produced 
by  shock  acceleration  ahead  of  fast  coronal  mass  ejections 
(CMEs;  Reames  1999;  Cliver  2009).  However,  the  intensities 
and  timescales  of  the  SEP  events  show  large  variations  even 
within  limited  ranges  of  solar  longitude  source  regions  and  CME 
speeds  and  widths  (Gopalswamy  et  al.  2004;  Kahler  2005). 

The  SEP  event  profiles  observed  at  1  AU  depend  on  a  number 
of  unknown  and  often  unobservable  factors,  such  as  the  detailed 
interaction  between  the  CME-driven  shock  and  the  SEP  seed- 
particle  population  (Tylka  &  Lee  2006;  Klecker  et  al.  2007).  A 
relatively  unexplored  area  is  the  role  of  the  coronal  environment 
into  which  the  fast  CME  propagates.  Gopalswamy  et  al.  (2002, 
2004)  found  that  SEP  event  peak  intensities  were  substantially 
enhanced  when  the  shock-driving  CME  was  preceded  within 
24  hr  by  another  CME  of  width  >60°  from  the  same  or  a 
nearby  solar  source.  Their  result  was  confirmed  by  Kahler  & 
Vourlidas  (2005),  who  also  found  that  CMEs  producing  more 
intense  SEP  events  were  more  massive  and  more  likely  to  be 
streamer-blowout  CMEs.  The  reasons  for  these  differences  are 
not  clear. 

Observations  of  SEP  events  within  various  solar  wind  (SW) 
structures  at  1  AU  give  little  indication  of  the  importance  of  those 
structures  for  SEP  events.  The  frequency  and  sizes  of  SEP  events 
(Kahler  2004)  as  well  as  their  timescales  (Kahler  2005)  in  fast 
SW  streams  are  comparable  to  those  in  slow  SW  streams.  More 


recently,  the  SEP  elemental  composition  observed  at  5-10  MeV 
nucleon-1  has  been  found  not  to  differ  between  the  two  kinds 
of  streams  (Kahler  et  al.  2009).  SEPs  can  also  propagate  to  very 
high  (>70°)  latitudes  (Sanderson  2004;  Malandraki  et  al.  2009) 
and  result  in  impressively  spatially  uniform  distributions  across 
broad  (>160°)  longitude  ranges  during  event  decay  phases 
(McKibben  et  al.  2003;  Reames  et  al.  1997). 

Since  the  SW  high-speed  streams  originate  from  coronal 
holes  (CHs;  Wang  2009),  these  results  suggest  that  CHs  are 
not  a  significant  factor  in  the  properties  of  SEP  events  observed 
at  1  AU.  Contrary  to  this  expectation,  however,  Kunches  & 
Zwickl  (1999)  found  that  the  times  from  flare  X-ray  maxima  to 
1  AU  onsets  of  >25  MeV  proton  events  during  1976-1997  were 
typically  longer  by  periods  of  hours  when  a  CH  lay  between  the 
CME  source  region  and  the  1  AU  magnetic  footpoint.  Their 
result  was  perhaps  inconsistent  with  later  studies  of  front-side, 
west-hemisphere  fast  halo  CMEs  by  Shen  et  al.  (2006,  2010). 
For  each  CME  they  looked  for  an  associated  1  AU  E  >  10  MeV 
or  >50  MeV  SEP  event  and  for  the  area,  shortest  distance,  and 
relative  location  of  any  large  CH  near  the  CME  source  region. 
They  found  that  the  CH  proximities,  areas,  and  locations  relative 
to  the  CME  source  and  Earth  magnetic  connection  were  not 
factors  in  the  occurrence  of  associated  SEP  events.  Shen  et  al. 
(2006,  2010)  considered  whether  SEP  events  occurred  for  given 
fast  halo  CMEs  while  Kunches  &  Zwickl  (1999)  compared 
onset  intervals  for  given  SEP  events,  so  their  results  are  not 
directly  comparable.  However,  the  two  works  suggest  opposite 
answers  to  the  basic  question  of  whether  CHs  play  any  role  in  the 
properties  of  SEP  events,  and  the  question  remains  controversial. 

How  can  a  CH,  or  its  high-speed  stream,  be  a  factor  in  the 
characteristics  of  an  SEP  event?  The  implication  in  the  Kunches 
&  Zwickl  (1999)  study  was  that  either  the  SEPs  must  propagate 
from  the  shock  source  across  the  intermediate  CH  field  lines  to 
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Figure  1.  Schematic  showing  how  the  pressure  of  magnetic  fields  might  deflect 
the  propagation  direction  of  a  CME.  The  CME  and  preceding  shock  (red  lines) 
that  would  have  propagated  toward  the  observer  (blue  dot)  are  deflected  eastward 
(blue  lines)  by  the  fields  extending  from  the  CH  lying  west  of  the  CME.  Adapted 
from  Gopalswamy  et  al.  (2010). 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

reach  those  connecting  to  Earth  or  the  speed  or  strength  of  the 
shock  producing  the  SEPs  is  diminished  as  it  crosses  the  field 
lines  extending  from  the  CH.  The  latter  possibility  of  diminished 
SEP  production  in  weaker  shocks  traversing  the  higher- Alfven 
speed  regions  characteristic  of  the  CH  (Kahler  &  Reames  2003) 
also  applies  to  coronal  field  regions  adjacent  to  the  CME. 

7.2.  CH  Deflections  ofCMEs 

Here  we  explore  an  alternative  possibility  that  the  CH  fields 
affect  coronal/interplanetary  SEP  production  by  latitudinally 
and/or  longitudinally  deflecting  the  CME  and  its  preceding 
shock  envelope.  This  work  is  motivated  by  the  recent  discovery 
(Gopalswamy  et  al.  2009)  that  some  disk-center  CMEs  produced 
shocks  at  1  AU  without  accompanying  in  situ  observations  of 
the  associated  interplanetary  CME  (ICME)  drivers.  Most  such 
driverless  shocks  occur  only  from  CMEs  near  the  solar  limbs, 
but  these  disk-center  CMEs  were  located  adjacent  to  CHs  which 
deflected  the  CMEs  away  from  the  Sun-Earth  line,  by  angles 
as  large  as  ~50°  (Gopalswamy  et  al.  2009,  2010).  The  CH 
deflection  effect  was  quantified  by  a  CH  influence  parameter 
(CHIP),  which  showed  statistically  that  disk-center  CMEs  pro¬ 
ducing  shocks  with  observed  drivers  were  negligibly  influenced 
by  CHs.  An  extended  statistical  comparison  (Mohamed  et  al. 
2012)  of  disk-center  CMEs  producing  the  six  driverless  shocks, 
shocks  with  observed  magnetic  clouds,  and  shocks  with  ICMEs 
found  the  CHIPs  to  increase  in  the  order  of  magnetic  cloud, 
ICME,  and  driverless  shocks.  This  result  suggested  magnetic 
flux  ropes  at  the  origin  of  all  the  disk-center  CMEs,  observable 
only  when  CH  deflections  are  small. 

Recent  observations  (Shen  et  al.  2011;  Gui  et  al.  2011;  Liu 
et  al.  2010;  Byrne  et  al.  2010;  Lugaz  et  al.  2011)  show  that 
these  CME  deflections  must  happen  within  several  RQ ,  where 
the  shock  is  expected  to  be  strongest.  The  magnetic  connectivity 
of  the  SEP  acceleration  region  to  the  observer  may  change  from 
a  favorable  to  an  unfavorable  one  or  vice  versa,  as  indicated  in 
the  schematic  of  Figure  1 .  Displacements  of  SEP  acceleration 
regions  suggest  a  new  SEP  environmental  factor  worthy  of 
further  exploration.  Latitudinal  deflections  of  CMEs  up  to  ~  40°, 
prominently  toward  the  equator  around  solar  minimum,  have 
been  studied  with  Solar  and  Heliospheric  Observatory /Large 


Angle  and  Spectrometric  Coronagraph  Experiment  (SOHO/ 
LASCO;  Gopalswamy  &  Thompson  2000;  Cremades  et  al. 
2006;  Wang  et  al.  201 1),  STEREO/ SECCHI  (Kilpua  et  al.  2009; 
Byrne  et  al.  2010;  Shen  et  al.  2011;  Zuccarello  et  al.  2012),  and 
SMEI  (Kahler  &  Webb  2010)  observations.  The  interpretation 
is  that  around  solar  minimum  the  CMEs  are  guided  equatorward 
by  strong  polar  CH  magnetic  fields.  This  effect  is  also  obvious 
in  the  rolling  motions  of  eruptive  prominences  erupting  in  the 
vicinity  of  CHs  (Panasenco  et  al.  2011,  2012).  Longitudinal 
CME  deflections  are  also  observed.  A  westward  motion  during 
the  acceleration  phase  found  in  SECCHI  CMEs  by  Liu  et  al. 
(2010)  was  attributed  to  coronal  corotation  of  CME  fields, 
although  adjacent  high-speed  streams  were  invoked  earlier  by 
Wang  et  al.  (2004)  to  explain  an  inferred  systematic  eastward 
deflection  of  fast  CMEs.  However,  the  physical  reality  of  CME 
deflections  by  tens  of  degrees  in  both  latitude  and  longitude  is 
now  well  established  by  both  models  and  observations. 

Modeling  of  the  magnetic-field  environment  has  shown  that 
CME  deflections  are  consistent  with  the  directions  and  magni¬ 
tudes  of  the  magnetic-field  energy-density  gradients  (Shen  et  al. 
2011;  Gui  et  al.  2011).  Furthermore,  slow  CMEs  are  deflected 
more  than  fast  CMEs  (Wang  et  al.  2004;  Gui  et  al.  2011).  An 
early  version  of  a  CH  magnetic  force  vector  F ,  derived  from  CH 
area  and  distance  from  the  CME  source  and  directed  against 
CMEs,  was  consistent  with  observed  latitudinal  CME  deflec¬ 
tions  (Cremades  et  al.  2006).  Gopalswamy  et  al.  (2009)  mod¬ 
ified  F  to  include  the  average  CH  magnetic  field  B  as  a  linear 
factor  of  their  CHIP  (Mohamed  et  al.  2012).  Gopalswamy  et  al. 
(2010)  then  redefined  F  to  a  pressure  term  with  a  B 2  factor, 
finding  that  F  was  larger  for  the  CMEs  of  driverless  shocks  than 
for  the  shocks  with  detected  drivers  at  1  AU. 

The  relevance  of  CME  deflections  to  SEP  events  rests  on 
two  considerations.  The  shock  flank  on  the  CH  side  may  be 
weakened  because  CHs  generally  have  higher- Alfven  speeds 
as  discussed  above.  The  second  is  the  question  of  shock 
connectivity  to  field  lines  well  connected  to  1  AU  observers. 
A  poorly  connected  (say  from  a  site  near  the  disk  center)  CME 
might  become  well  connected  if  the  CME  nose  is  deflected  to  the 
west  by  a  CH  located  between  the  disk  center  and  the  eruption 
site.  On  the  other  hand,  the  connectivity  may  be  diminished  if  the 
CH  is  located  to  the  west  of  the  eruption  region,  as  in  Figure  1 . 
The  presence  of  one  or  more  CHs  near  the  CME  source  region 
could  modify  the  SEP  production  in  either  case,  and  in  this  work 
we  seek  evidence  for  this  effect. 

2.  DATA  ANALYSIS 
2.7.  Selection  of  SEP  and  CH  Data 

For  the  analysis  here,  we  use  timing  parameters  and 
peak  intensities  of  20  MeV  SEP  events  measured  with  the 
EPACT  instrument  on  the  Wind  spacecraft  and  associated  with 
SOHO /CASCO  CMEs.  The  event  list  of  Table  1  of  Kahler 
(2005)  has  been  updated  through  2008.  Since  any  effects  of 
CHs  on  CME  and  shock  propagation  and  associated  SEP  ac¬ 
celeration  should  be  most  apparent  during  the  rise  phase  of  the 
SEP  event,  we  select  for  analysis,  besides  the  peak  SEP  inten¬ 
sity  Ip,  the  SEP  onset  time  T0 ,  defined  as  the  time  from  CME 
onset  to  SEP  onset  at  Wind ,  and  TR,  the  time  from  SEP  onset 
to  half-maximum  of  the  SEP  peak  intensity  (Figure  2).  These 
SEP  timescales  can  range  up  to  an  order  of  magnitude  (Table  2, 
Kahler  2005),  even  within  relatively  narrow  source  longitude 
ranges.  Here  we  look  for  any  effects  of  CHs  on  the  SEP  event 
timescales  and  intensities. 
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Figure  2.  Plot  of  2  and  20  MeV  proton  intensities  from  the  Wind /EP ACT 
instrument  during  the  SEP  event  of  2002  July  16.  The  horizontal  bars  show  To, 
the  time  from  CME  onset  to  SEP  onset  at  Wind ;  Tr,  the  rise  time  from  onset  to 
half  the  peak  intensity;  and  To,  the  event  duration  to  half  peak  or  shock.  Only 
Tq  and  Tr  are  used  in  the  analysis  here. 


Table  1 

The  CH  F  Parameters  for  the  2004  November  7  SEP  Event 


Location 

Areaa 

Ave  B 

Distance13 

F 

FPAC 

Flat 

Flong 

S49W06 

1.3e+ll 

2.1 

6.8e+05 

5.9e-01 

354° 

5.6e-02 

5.9e-01 

S05E11 

7.9e+09 

3.6 

3.9e+05 

1.9e-01 

298° 

1.7e-01 

9.3e-02 

N37E36 

3.6e+10 

-4.9 

6.4e+05 

4.2e-01 

231° 

3.3e-01 

— 2.6e-01 

N38E61 

2.0e+10 

-10.0 

8.4e+05 

2.9e-01 

233° 

2.3e-01 

—  1.7e-01 

N63E42 

5.8e+10 

-4.1 

7.7e+05 

4.0e-01 

205° 

1.7e-01 

— 3.6e-01 

Total  F 

9.6e-01 

264° 

9.6e-01 

—  1.0e-01 

Notes. 

a  Area  in  units  of  km2. 

b  Distance  in  units  of  km. 

c  F  position  angle  measured  ccw  from  north. 

CHs  and  CME  source  regions  are  best  observed  at  solar  cen¬ 
tral  meridian  (CM),  so  we  select  a  source  longitude  range  of 
E20°  to  W20°,  which  yields  41  SEP  events  near  CM  with 
relatively  well-observed  CH  boundaries  and  intensities.  We 
identify  the  CME  source  regions  by  the  associated  flare  ac¬ 
tive  region  locations,  and  use  the  SOHO  EIT  284  A  images  to 
define  the  CH  boundaries  and  the  SOHO  MDI  images  to  deter¬ 
mine  the  average  CH  magnetic  fields.  Following  Gopalswamy 
et  al.  (2010),  we  calculate  for  each  CME  the  influence  parameter 
vector  F,  defined  as 

F  =  B2A/r2, 

where  A  is  the  area  of  the  CH,  B  is  the  average  magnetic  field  in 
the  CH,  and  r  is  the  CH  distance.  F  for  each  CH  is  directed  from 
the  CH  center  to  the  CME  source  region,  and  the  total  F,  in  units 
of  Gauss  squared,  is  the  vector  sum  of  F  over  all  individual  CHs. 
The  position  angle  FPA  is  the  assumed  direction  in  which  the 
CH  fields  are  deflecting  the  CME,  as  indicated  for  an  example 
SEP  event  in  Figure  3.  The  parameters  for  each  of  the  five  CHs 
and  their  F  vectors  in  this  case  are  given  in  Table  1 ,  where  the  last 
two  columns  give  the  latitudinal  Flat  (positive  is  northward)  and 
longitudinal  Flong  (positive  is  westward)  vector  components. 
The  tail  of  the  arrow  showing  F  is  placed  on  the  CME  source 
region  in  Figure  3  to  indicate  the  expected  deflection  direction 
of  the  CME. 
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Figure  3.  Full-disk  EIT  284  A  image  on  2004  November  7  showing  boundaries 
of  the  five  CHs  used  to  calculate  the  component  F  vectors  of  Table  1 .  The  red 
dotted  lines  show  the  connecting  distances  r  and  directions  between  the  CHs  and 
the  CME  source  region  for  each  component.  The  position  angle  of  the  summed 
F  vector  (red  arrow)  is  264°  and  is  placed  on  the  active  region  source  of  the 
CME  at  N09W17  to  show  the  expected  CME  deflection  direction. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

Since  the  presumed  influence  of  CHs  on  SEP  events  is 
produced  by  deflections  of  CMEs,  we  also  look  at  the  CME 
propagation  directions  given  in  the  CDAW  LASCO  CME 
catalog.  Of  the  41  CMEs,  30  are  full  halos  spanning  360° 
of  width,  but  each  CME  speed  is  determined  from  a  leading- 
edge  direction,  the  measurement  position  angle  (MPA),  selected 
from  the  first  CME  images.  These  directions  are  usually,  but 
not  always,  near  the  central  angle  of  the  CME  span.  In  our 
examination  of  the  movies  of  subtracted  C2  LASCO  images, 
we  found  six  cases  of  significant  disagreements  ranging  from 
25°  to  105°  between  the  catalog  MPA  and  our  determination  of 
the  primary  CME  PA.  We  used  our  6  PAs  and  the  catalog  MPA 
values  for  the  other  35  CMEs  in  the  data  analyses. 

2.2.  The  Statistical  Comparisons 

2.2.1.  The  CH  Deflection  Correlations 

In  this  statistical  study,  we  look  for  correlations  between  pairs 
of  SEP  event  and  associated  CH  F parameters.  The  expectation  is 
that  CME  deflections  either  westward  or  equatorward  from  high- 
latitude  source  regions  will  result  in  more  favorable  magnetic 
connections  of  shock  accelerated  SEPs  to  the  field-line  con¬ 
nections  to  the  Earth,  resulting  in  shorter  SEP  event  timescales 
and/or  higher  peak  intensities.  For  this  purpose,  we  deal  sepa¬ 
rately  with  the  Flong  and  Flat  components  of  F.  We  redefine  Flat 
to  be  positive  when  pointing  from  the  solar  source  region  toward 
the  ecliptic  and  negative  when  pointing  poleward.  The  sign  of 
Flat,  i.e.,  whether  the  CME  was  being  deflected  equatorward  or 
poleward,  was  determined  after  first  calculating  the  solar  source 
latitude  relative  to  the  ecliptic  plane  using  the  solar  B  angle. 
There  were  27  westward  and  14  eastward  Flong  components 
and  16  ecliptic  ward  and  25  poleward  Flat  components  among 
the  41  total  events. 

We  do  a  series  of  comparisons  to  look  for  any  significant 
correlations  between  the  SEP  event  and  CH  F  parameters.  The 
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Figure  4.  Logs  of  the  SEP  event  onset  times  To  vs.  the  associated  CME  MPA 
directions.  Most  MPAs  are  in  the  western  hemisphere  (180°  <  MPA  <  360°), 
and  there  is  a  slight  but  not  significant  negative  correlation  of  To  with  MPA 
(Table  2),  suggesting  a  marginally  shorter  To  for  westward  propagating  and 
better  connected  CMEs. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


Table  2 

Correlation  Coefficients  between  the  SEP  Events  and  CH  or  CME  Parameters 


SEP/CME  (No.  of  Events) 

Log  Ipeak 

Log  T0 

Log  Tr 

90%  CC 

Log  F  (41) 

0.24 

0.17 

0.11 

0.26 

Flong  (41) 

0.00 

0.14 

0.06 

0.26 

Flat  (41) 

0.09 

-0.21 

-0.06 

0.26 

Log  West  Flong  (27) 

-0.07 

0.25 

-0.10 

0.32 

Log  East  Flong  (14) 

0.10 

0.22 

-0.00 

0.46 

Log  Equat  Flat  (16) 

0.27 

-0.28 

0.05 

0.43 

Log  Pole  Flat  (25) 

-0.03 

-0.35 

-0.22 

0.34 

MPA  (41) 

0.06 

-0.21 

-0.17 

0.26 

correlation  coefficients  (CC)  are  presented  in  Table  2  along  with 
the  CC  required  for  a  modest  90%  probability  of  correlation 
(Bevington  &  Robinson  2003)  in  the  last  column.  If  we  consider 
that  the  directions  rather  than  the  magnitude  of  F  would  be 
important  for  ordering  SEP  events  by  producing  better  or  worse 
magnetic  connections  of  SEP  sources,  then  any  correlations  of 
SEP  parameters  with  log  F,  given  in  the  first  line  of  Table  2, 
would  not  be  expected.  We  compared  the  SEP  parameters  with 
linear  values  of  Flong  and  Flat,  which  ranged  over  negative  and 
positive  values,  and  then  separately  with  logs  of  the  east,  west, 
equatorward,  and  poleward  components,  respectively  in  rows 
four  to  seven  of  Table  2.  The  statistics  are  now  diminishing, 
but  the  only  hints  of  correlations  seem  to  go  in  an  unexpected 
direction,  such  as  that  (CC  =  —0.35,  90%  significance  for  25 
events)  of  the  shorter  onset  times  To  with  increasing  strengths 
of  the  poleward  deflections  of  Flat  or  with  decreasing  strengths 
of  the  westward  deflections.  The  primary  result  here  is  that  we 
find  no  significant  correlations  between  the  SEP  event  and  CH 
F  parameters,  as  shown  in  Table  2. 

2.2.2.  The  SEP  Parameters  versus  the  CME  MPAs 

The  analysis  of  the  previous  section  was  based  on  the  two 
assumptions  that  the  SEP  event  parameters  might  be  dependent 
on  the  propagation  direction  of  the  CME  and  that  the  CME 
propagation  would  be  deflected  by  the  CH  magnetic  fields.  The 
lack  of  any  correlation  between  the  CH  F  vectors  and  the  SEP 
event  parameters  now  raises  two  questions.  First,  do  the  SEP 
parameters  depend  on  the  CME  propagation  directions,  and 
second,  do  the  CH  fields  significantly  deflect  the  CMEs?  We 
can  use  our  data  set  of  41  SEP  events  and  CMEs  to  investigate 
both  questions. 
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Figure  5.  Source  flare  PA  (fPA)  correlates  with  the  observed  CME  MPA  at 
CC  =  0.47  for  the  36  events.  A  modified  plot  (see  the  text)  yields  a  CC  =  0.91, 
suggesting  a  radial  propagation  of  the  CME  from  the  flare  site. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


We  might  expect  that  the  SEP  timescales  are  shortest  and  peak 
intensities  highest  for  westward  CME  propagation  directions 
and  reversed  for  unfavorable  eastward  CME  directions,  resulting 
in  extreme  values  at  CME  MPAs  centered  on  ~270°  and 
~90°,  respectively.  However,  plotting  each  of  the  SEP  event 
parameters  against  the  CME  MPAs,  we  find  no  significant 
correlations,  as  shown  on  the  last  line  of  Table  2.  While 
the  plot  of  log  Tq  against  the  CME  MPA  in  Figure  4  is 
slightly  suggestive  of  the  expected  shorter  values  at  ~270°, 
the  characteristically  broad  scatter  of  To  across  the  MPA 
values  indicates  an  independence  of  the  SEP  parameters  on 
the  directions  of  CME  propagation  from  CM  source  regions. 


2.2.3.  Significance  of  the  Flare  Locations 
and  F  Vectors  for  CME  MPAs 

The  next  question  is  how  effectively  the  CH  magnetic  force  F 
moderates  the  CME  MPAs.  We  first  converted  the  heliographic 
latitude  coordinates  of  the  associated  flare  locations  to  ecliptic 
latitude  coordinates  by  correcting  for  the  solar  B  angle  and  then 
calculated  a  solar  PA  relative  to  the  sub-Earth  point  for  each  of 
the  associated  36  flare  sources  lying  outside  a  15°  radius  of  the 
sub-Earth  point.  The  source  flare  PA  (fPA)  correlates  with  the 
observed  MPA  at  CC  =  0.47,  well  above  the  99%  significant 
level,  as  shown  in  Figure  5.  However,  if  we  are  mindful  of  the 
discontinuity  of  adjacent  PAs  from  0°  to  360°  and  reorganize  the 
plot  by  converting  two  fPA  values  in  the  upper  left  corner  from 
5°  and  21°  to  365°  and  381°,  respectively,  and  a  third  value  in 
the  lower  right  of  the  plot  from  356°  to  —4°,  then  the  CC  rises 
to  0.91,  well  exceeding  the  99.9%  significance  level.  Thus  the 
fPA  is  a  reasonable  guide  to  the  MPA  of  the  CME  propagation, 
which  suggests  a  generally  radial  propagation  of  the  CMEs  from 
their  source  regions  independent  of  any  CH  effects. 

Despite  the  excellent  agreement  of  the  CME  MPAs  with  the 
fPAs  in  Figure  5,  we  can  still  compare  the  PAs  of  the  F  vectors 
(FPAs)  with  the  41  observed  MPAs  to  look  for  a  correlation 
that  would  indicate  CME  deflections  in  the  directions  of  the  F 
vectors.  The  plot  of  Figure  6  has  a  correlation  of  CC  =  0.36, 
significant  at  the  98%  significance  level.  There  is  more  scatter 
in  this  plot  than  that  of  the  MPA  versus  fPA  of  Figure  5,  which 
suggests  that  the  CH  magnetic  fields  may  play  some  role  in 
deflecting  CME  propagation  directions. 
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MPA  (deg) 

Figure  6.  Comparison  of  the  41  CH  FPA  values  with  the  observed  MPAs  also 
shows  a  significant  correlation  consistent  with  an  effect  of  the  CH  fields  on  the 
direction  of  CME  propagation. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 
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Figure  7.  DPA  may  not  be  a  useful  parameter  to  judge  the  effect  of  CH  fields  on 
CMEs,  as  illustrated  in  the  2001  April  9  event  in  the  same  format  as  Figure  3. 
The  flare  source  site  15°  south  of  the  sub-Earth  point  implies  a  southward  CME 
propagation,  which  may  have  been  modified  by  the  southern  polar  CH  toward  a 
lower  latitude  course,  but  not  sufficient  to  change  the  MPA  of  21 1°. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

2.2.4.  The  DPA  as  a  Further  Test  of  CH  Field  Effectiveness 

Gopalswamy  et  al.  (2009)  calculated  the  difference  between 
the  MPA  and  the  FPA,  defined  as  the  DPA  (At/r  in  Mohamed 
et  al.  2012),  a  measure  of  how  well  the  CME  propagation  angle 
correlated  with  that  expected  from  an  effective  CH  F  vector. 
In  their  examination  of  shocks  without  observed  driver  ICMEs, 
they  found  a  better  agreement,  i.e.,  smaller  DPA  values,  for 
those  CM  CMEs  not  observed  at  1  AU  than  for  those  that  were 
observed.  This  result  was  consistent  with  the  idea  that  CMEs 
not  observed  at  1  AU  were  effectively  deflected  by  the  CH 
magnetic  forces.  However,  the  DPA  may  not  always  be  a  useful 
parameter  to  judge  effects  of  CH  fields  on  CMEs,  as  illustrated 
for  the  2001  April  9  event  shown  in  the  EIT  image  of  Figure  7. 
The  DPA  for  that  event  is  144°,  a  large  value  resulting  from  an 
MPA  of  211°  and  an  FPA  of  355°.  A  second  SEP  event  on  the 
following  day  with  the  same  source  region  had  a  similarly  high 
DPA  of  156° .  Since  the  CME  source  region  lies  south  of  the  sub- 


Figure  8.  DPA  vs.  the  CME  speeds  for  the  41  SEP  events.  Larger  DPA  values 
indicate  worse  agreement  between  the  CME  propagation  direction  MPA  and  the 
CH  force  vector  F  direction.  This  result,  significant  at  the  90%  level,  supports 
the  concept  the  CH  deflection  force  is  less  effective  for  faster  CMEs,  presumed 
to  have  larger  momenta  than  slower  CMEs. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

Earth  point,  we  expect  that  the  MPA  should  be  southward,  near 
180°,  as  observed.  The  expected  effect  of  the  magnetic  fields 
from  the  dominant  southern  CH  should  be  to  deflect  the  CME 
equatorward,  but  if  that  deflection  occurred,  it  was  not  sufficient 
to  change  the  resulting  southward  projection  of  the  CME. 

With  this  caveat  of  a  projection  effect  in  mind,  we  have 
compared  DPA  =  |MPA  -  FPA|  with  log  F  to  see  whether, 
as  the  CH  magnetic  force  increases,  the  DPA  becomes  smaller, 
indicating  better  agreement  with  the  expected  deflection  PA  and 
resulting  in  a  negative  correlation  between  the  two  parameters. 
However,  we  instead  calculate  a  small  positive  CC  =  0. 15,  below 
the  significance  value  of  0.26,  but  opposite  to  the  expected  result 
for  significant  CH  field  effects  on  CMEs.  A  test  for  the  projection 
effect  is  to  look  for  a  trend  for  the  CME  sources  with  angular 
separations  farthest  from  the  sub-Earth  point  to  have  the  largest 
DPAs,  as  we  discussed  for  the  example  of  Figure  7.  However, 
we  find  CC  =  0.02  for  DPA  versus  source  separation,  consistent 
with  no  effect  of  the  CH  F  on  the  CME  MPA. 

Another  approach  is  to  ask  whether  the  faster  CMEs  are  less 
deflected  than  slower  CMEs,  with  presumed  smaller  momenta, 
by  the  CH  fields,  as  suggested  by  Wang  et  al.  (2004)  and  Gui 
et  al.  (201 1).  For  this  purpose  we  compared  the  DPA  with  CME 
speed  and  found  a  CC  =  0.28,  a  barely  significant  value  at  the 
90%  significance  level,  but  clearly  supportive  of  an  effective 
CH  force  on  CME  propagation  (Figure  8).  Thus,  the  DPA 
correlation  analyses  give  conflicting  results,  increasing  with  F, 
but  decreasing  with  CME  speed.  If  we  give  more  weight  to  the 
latter  result,  then  the  DPA  comparisons  provide  weak  evidence 
for  significant  deflections  of  CMEs  by  the  CH  fields. 

3.  SUMMARY  AND  DISCUSSION 
3.1.  Effects  of  CH  Fields  and  CME  Propagation  on  SEP  Events 

The  recent  result  that  CH  magnetic  fields  can  significantly 
deflect  shock-producing  CMEs  from  near  CM  (Gopalswamy 
et  al.  2009,  2010)  raised  the  possibility  that  these  deflections 
may  have  significant  consequences  for  the  SEP  events  produced 
by  particle  acceleration  in  the  shocks  driven  by  CMEs.  We  used 
a  sample  of  41  20  MeV  SEP  events  with  associated  CMEs  from 
CM  locations  to  look  statistically  for  systematic  variations  in 
SEP  onset  T0  and  rise  TR  times  and  peak  intensities  Ip  resulting 
from  CME  deflections.  An  ad  hoc  force  parameter  F connects  the 
CH  properties  of  magnetic-field  intensity,  area,  and  proximity  to 
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the  expected  CME  deflection.  The  question  is  whether  the  broad 
range  of  SEP  event  peak  intensities  (Kahler  2001)  and  timescales 
(Kahler  2005)  can  be  organized  in  any  way  by  the  use  of  the 
calculated  F  vector.  We  first  sought  a  direct  connection  between 
SEP  parameters  and  the  magnitude  or  direction  of  the  F  vectors, 
which  were  resolved  into  equatorward/poleward  and  east/ west 
components.  We  found  no  statistically  significant  correlations 
between  these  two  sets  of  parameters,  as  indicated  in  Table  2. 

The  independence  of  SEP  event  parameters  from  expected 
CME  deflections  led  us  to  ask  whether  the  SEP  parameters 
depend  on  the  CME  propagation  directions  as  defined  by  their 
MPAs.  SEP  event  onset  and  rise  times  are  dependent  on  the 
longitudes  of  the  solar  source  regions  (Kahler  2005),  but  with 
a  selection  of  source  regions  from  a  longitudinal  band  only  40° 
wide,  we  might  expect  that  their  associated  CME  propagation 
directions  would  provide  some  of  the  longitudinal  organization 
lost  by  confining  our  source  regions  to  the  40°  wide  band.  That 
this  is  not  the  case  is  obvious  from  Figure  4  in  which  we  find 
no  significant  difference  between  eastward  (MPA  ~  90°)  and 
westward  (MPA  ~  270°)  CMEs  for  T0.  The  implication  here 
is  that  even  if  CH  fields  are  deflecting  CMEs,  those  deflections 
are  not  relevant  to  the  SEP  events  produced  by  the  CME-driven 
shocks. 

3.2.  CME  Propagation  Directions  and  CH  Fields 

We  also  asked  whether  we  find  evidence,  independent  of  the 
SEP  events,  for  an  impact  of  the  CH  magnetic  fields  in  deflect¬ 
ing  CMEs.  Comparing  first  the  CME  propagation  directions 
expected  only  on  the  basis  of  the  associated  fPAs,  we  found 
the  good  agreement  shown  in  Figure  5.  While  this  result  could 
suggest  an  essentially  undeflected  radial  propagation  for  those 
CM  CMEs,  we  also  found  a  modest  but  significant  correlation 
between  the  FPA  and  MPA  (Figure  6),  consistent  with  some  CH 
effect  in  deflecting  the  CME  propagation  direction.  From  their 
proximity  to  the  Sun  center,  we  can  expect  substantial  scatter  in 
the  MPAs  of  the  CMEs  of  this  study,  which  may  contribute  to 
the  scatter  of  Figure  6.  The  MPA  provides  us  with  only  an  az¬ 
imuthal  direction  projected  on  the  plane  of  the  sky,  so  significant 
CME  deflections  of  propagation  from  radial  can  go  undetected, 
as  illustrated  in  Figure  7  showing  an  FPA  almost  opposite  in  di¬ 
rection  from  the  observed  MPA.  This  could  explain  why  we  find 
a  weak  increase  of  DPA  with  increasing  magnitude  of  F,  a  result 
contrary  to  the  expectation  that  stronger  CH  fields  F  should  be 
more  effective  in  deflecting  CMEs.  STEREO/ SECCHI  obser¬ 
vations  (e.g.,  Liu  et  al.  2010)  will  be  required  to  determine  the 
actual  deflections  from  radial  propagation  that  might  result  from 
CH  fields. 

An  alternative  and  complementary  approach  to  look  for  evi¬ 
dence  of  CME  deflections  is  to  compare  latitudinal  deviations 
of  limb  CMEs  with  their  calculated  CH  latitudinal  F  vectors. 
This  approach  has  obvious  limitations  in  determining  appropri¬ 
ate  CH  properties,  but  avoids  the  projection  limitations  of  this 
study  indicated  in  Figure  7.  The  results  of  such  a  study,  shown  in 
Figure  6  of  Cremades  et  al.  (2006),  indicate  a  good  correlation 
between  CME  deviations  and  F  vectors. 

As  a  second  test  of  the  effectiveness  of  CH  fields  in  deflecting 
CMEs,  we  followed  the  assertions  of  Wang  et  al.  (2004)  and 
Gui  et  al.  (2011)  that  slower  CMEs  should  be  more  susceptible 
than  faster  CMEs  to  the  deflecting  forces  of  CH  fields.  Our 
comparison  of  DPA,  the  angular  difference  between  MPA  and 
FPA,  with  CME  speed  shows  good  agreement  with  that  idea. 
The  results  shown  in  Figures  6  and  8  further  support  the 
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concept  of  CME  deflections  by  CH  fields  (Cremades  et  al.  2006; 
Gopalswamy  et  al.  2009). 

3.3.  Implications  and  Limitations  of  the  Study 

A  fraction  of  the  variability  of  SEP  event  timescales  is 
dependent  on  solar  source  longitude  (Kahler  2005),  but  we 
find  no  evidence  that  within  a  confined  longitude  range  SEP 
event  parameters  depend  on  the  CME  propagation  directions 
defined  by  their  MPAs.  We  had  expected  from  Figure  1  that 
CME  propagation  directions  should  be  at  least  a  contributing 
factor  for  the  access  of  SEPs  to  Earth-connected  interplanetary 
field  lines  and  would  be  manifested  in  the  statistical  distributions 
of  the  SEP  event  timescales  and  perhaps  peak  intensities.  In  an 
earlier  study  Kahler  et  al.  (2000)  found  that  SEP  peak  intensities 
were  independent  of  CME  latitudes  for  a  set  of  17  SEP  events 
limited  to  an  SMMC/P  CME  speed  range  of  650-850  km  s-1 
and  source  longitudes  of  >W20°.  Our  work  with  CM  sources 
allows  us  to  examine  a  full  range  of  CME  MPAs  and  extend  the 
comparisons  to  include  SEP  event  timescales.  The  lack  of  any 
SEP  event  dependence  on  the  MPAs  of  these  CMEs  is  consistent 
with  SEP  production  in  CME-driven  shocks  much  broader  than 
those  suggested  in  the  schematic  of  Figure  1 .  SEP  events  (Kahler 
&  Reames  2003)  and  type  II  radio  shock  bursts  (Gopalswamy 
et  al.  2001)  are  associated  only  with  CME  widths  >60°,  but 
the  inferred  angular  extents  of  shocks  at  1  AU  can  reach  >100° 
(Gopalswamy  et  al.  2001 ;  de  Lucas  et  al.  201 1),  and  SEP  events 
have  been  observed  ~180°  from  their  source  regions  (Cliver 
et  al.  2005). 

Figures  6  and  8  show  correlations  supporting  the  idea  of  CH 
deflections  of  CMEs,  but  one  possible  bias  is  that  CHs  are  often 
located  between  the  CME  source  active  regions  and  Sun  center 
such  that  the  F  vectors  coincidently  align  with  the  CME  MPAs 
determined  primarily  by  the  CME  source  locations.  This  could 
then  give  rise  to  the  apparent  correlation  of  FPAs  with  MPAs  of 
Figure  6.  We  think  this  is  unlikely  since  F  is  typically  calculated 
from  multiple  CHs.  In  using  the  associated  flare  locations  as 
the  CME  sources  we  have  ignored  offsets  of  CME  sources 
from  active  region  centers  (Wang  &  Zhang  2007)  and  CME 
deflections  from  nearby  sunspot  fields  (Sterling  et  al.  2011). 
Our  study  was  limited  to  SEP  events  with  sources  within  20°  of 
CM  and  was  based  on  an  ad  hoc  force  vector  F  calculated  from 
the  properties  of  all  CHs  observed  on  the  Sun  at  the  time  of  each 
SEP  event.  We  can  not  rule  out  the  possibility  that  SEP  events 
with  origins  at  other  longitude  regions  may  in  fact  show  some 
CH  dependence,  but  our  result  implies  that  effects  of  CHs  are  not 
of  importance  for  forecasting  properties  of  SEP  events.  Finally, 
we  point  out  that  this  is  only  a  coronal  study,  which  does  not 
include  a  consideration  of  the  extension  of  the  CH  environment 
to  the  fast  SW  stream  regions  at  1  AU.  Effects  of  CHs  may  still  be 
important  for  SEPs,  but  only  in  the  interplanetary  medium  where 
the  magnetic  connectivity  of  the  observer  to  particular  solar 
source  regions,  and  CHs  in  particular,  may  play  a  significant 
role.  A  complementary  study  to  compare  these  SEP  events  with 
interplanetary  fields  and  their  solar  footpoints  will  be  done. 
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